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1.  INTRODUCTION 


It  is  very  much  possible  to  reshape  the  viscoelastodynamic  equations  into  a  form  that  resembles  or 
parallels  Maxwell’s  equation  to  achieve  an  analogy  or  a  complete  mathematical  equivalence;  hence,  a 
mathematical  parallelism  can  be  generated  between  electromagnetics  and  vibratory/acoustic  energy  waves. 
This  work  will  show  that  the  wave  equation  in  a  Maxwell  anisotropic-viscoelastic  solid  is  analogous  to 
the  two-dimensional  (2-D)  Maxwell  equations  describing  propagation  of  the  transverse  electromagnetic 
mode  in  anisotropic  media.  This  type  of  analogy,  i.e.,  between  the  process  of  conduction-static  induction 
through  dielectrics  and  viscosity  -  elasticity,  was  probably  assumed  initially  by  Maxwell  himself  (Eveiitt 
1975). 

This  phenomenological  analogy  can  be  utilized  in  many  ways.  Many  viscoelastodynamic  modeling 
codes  can  be  changed  to  simulate  electromagnetic  wave  propagation.  In  addition,  many  sets  of  solutions 
of  the  viscoelastic  problem  can  be  used  to  test  electromagnetic  codes.  Furthermore,  plane  propagation 
wave  theory  of  harmonic  waves  especially  involving  anisotropic-viscoelastic  media  apply  analogously  to 
electromagnetic  anisotropic  wave  propagation.  It  is  well  known  that  velocity  and  attenuation  anisotropy 
of  vibratory  waves  are  important  in  frequency-amplitude  vs.  time  analyses  conducted  on  the  Ml  Abrams 
tank  gas  turbine  engine  (Koijack  1995). 

2.  MAXWELL’S  EQUATIONS 

The  Maxwell  equations  can  be  written  as  (Chew  1990  and  Jordan  1950): 


VxE  =  - 


dB 

3t 


-M, 


(1) 


VxH  =  .^+J,  (2) 

where  E,  B,  H,  and  D  are  the  electric  field  intensity,  the  magnetic  flux  density,  the  magnetic  field  intensity 
and  the  electric  flux  density,  respectively,  and  J  and  M  are  the  electric  and  magnetic  current  densities, 
respectively.  Equations  (1)  and  (2)  constitute  six  scalar  equations  with  12  scalar  unknowns,  since  M  is 


1 


assumed  to  be  given  and  J  is  a  known  function  of  the  electric  field  as  stated  explicitly  by  Equation  (5) 
below.  The  six  additional  scalar  equations  are  the  constitutive  relations,  which  for  isotropic  media  can 
be  written  as 


D  =  e  ’E, 

(3) 

B  =  p  •  H, 

(4) 

where  e  and  p  are  the  permittivity  and  permeability  tensors,  respectively.  The  product  appearing  in 
Equations  (3)  and  (4)  is  simply  Cayley  multiplication.  Moreover,  the  current  density  is 


J  =  o  •  E  +  J,,  (5) 

where  c  is  the  conductivity  tensor  and  Jg  is  the  given  contribution  of  the  sources.  Substituting  the 
constitutive  relations  and  the  current  density  into  Equations  (1)  and  (2)  gives 

VxE  =  -p-.^+M,  (6) 

at 

VxH  =  a*E  +  e'.^+Jg.  (7) 


3.  THE  SONIC-WAVE  HELD  EQUATIONS 

The  elementary  equations  of  acoustics,  as  written  in  terms  of  particle  velocity  and  stress,  can  be 
expressed  in  terms  of  first  order  time  derivatives.  Cauchy’s  equations  can  be  written  as  (Auld  1990), 

V-T  =  p.^-F,  (8) 

at 


where 
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is  the  stress  tensor,  v  is  the  particle  velocity  tensor,  p  is  the  density,  F  is  the  body  force  tensor.  The  DEL 
(V)  operator  can  be  simply  expressed  for  rectangular  cartesian  coordinates  as. 


V  = 


d/dx 

0 

0 


0  0  0  d/dz  d/dy 

d/dy  0  d/dz  0  d/dx 

0  d/dz  d/dy  d/dx  0 


(10) 


The  strain  can  be  given  in  terms  of  the  displacement  u 


by, 


^  ~  t^xx’  ^y’  ^^z’ 


(11) 


such  that  e^x  =  dajdx,  e^y  =  (du^/dy  +  dUy/dx)/2,  etc.  Strain  and  particle  velocity  can  be  related  as. 


ttT  dS 
V*  'V  = 

A 


(12) 


The  sonic-electromagnetic  analogy  was  first  established  by  Auld  (1990)  by  using  a  3-D  Kelvin- Voigt 
Model, 


T 


c  dS 


(13) 


where  c^  and  %  are  the  (Kelvin-Voigt)  elasticity  and  viscosity  tensors,  respectively.  If  we  use  the  time 
derivative  of  the  strain  by  using  Equation  (12),  and  defining  the  tensor. 


it  can  be  deduced  that. 


at  9t 


(14) 


By  introducing  the  3-D  Maxwell  constitutive  relation  (Casula  and  Carcione  1992),  we  have, 


as  ^-1  ax 


-nM-T, 


(15) 


such  that  Cj^^  and  Tij^  are  the  elasticity  and  the  viscosity  tensors,  respectively.  This  relation  compares  very 
similarly  to  the  1-D  Maxwell  stress-strain  relation  (Ben-Menahem  and  Singh  1981).  If  we  drop  out  the 
strain  tensor  by  invoking  Equation  (12),  we  end  up  with  a  relation  similar  to  Greenfield  and  Wu  (1991) 
as. 


-1 


-v  =  11m 


T*Cm' 


at 


(16) 


If  a  compliance  tensor  can  be  defined  as 


along  with  the  complementary  compliance  tensor  of 


(17) 


(18) 


then  our  Equation  (16)  simply  reverts  to  the  expression. 


V'^-v  =  Xm  -T  +  s 


ar 


(19) 
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Generally,  the  analogy  does  not  mean  that  sonic-wave  and  electromagnetic  equations  illustrate  the 
exact  same  mathematical  problem.  In  actuality,  T  is  a  6-D  vector  and  E  is  a  3-D  vector.  In  addition, 
acoustics  involve  at  least  order  6  matrices  for  material  properties  and  electromagnetism  involves  merely 
order  3  matrices.  However,  a  version  of  complete  equivalence  can  be  adequately  established  in  the  2-D 
case  by  using  the  Maxwell  model. 

4.  SONIC-WAVE  AND  ELECTROMAGNETIC  ANALOGY 

By  considering  anisotropic  permittivity  and  conductivity  tensors  with  symmetry,  we  may  imagine  a 
medium  such  that 

Ell  0  £i3 

e  =  0  ^2  ®  (20) 

5i3  0 

and 

Oil  ®  ®13 

o=  0  022  0  .  (21) 

Oi3  0  033^ 

These  previously  mentioned  tensors  are  for  a  medium  having  the  vertical  axis  at  90°  to  the  symmetric 
plane.  It  is  well  known  that  a  coordinate  transformation  exists  that  fully  diagonalizes  these  symmetric 
tensors,  which  is  commonly  referred  to  as  the  principal  system  of  the  particular  medium  in  question; 
hence,  the  three  principal  components  of  the  tensors  ,  become  immediately  apparent  The  permeability 
tensor  is  usually  isotropic.  Here,  it  is  assumed  that  p  =  p  •  I  such  that  p  is  the  permeability  and  I  is  the 
identity  tensor  of  order  3. 

If  we  allow  the  material  properties  to  not  change  in  the  y-direction  and  only  allow  propagation  to 
occur  in  the  x-z  plane,  then  it  follows  that  E^,  H^^,  and  H^  are  the  only  components  under  consideration; 
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therefore,  if  there  are  no  electric  sources  within  the  media,  it  follows  that  the  field  equations  for  the 
transverse  electric  and  magnetic  (TEM)  field  take  on  the  form. 


an, 

ax 


aEj 

“ST 


(22) 


aH„ 

=  OiiEx+OigEj  +  eji 


+ 


(23) 


aHy 

ax 


aE^  aE^ 

®13  +  ®33  Ez  +  £i3 -gp  +  %3 • 


(24) 


According  to  Greenfield  and  Wu  (1991),  Equations  (22)  through  (24)  conceptualize  an  isotropic  model. 
In  addition,  if  we  have  uniform  properties  in  the  y-direction  as  in  the  field  of  sonics,  then  we  have  the 
situation  that  one  of  the  shear  waves  h^  its  own  differential  equation,  which  is  decoupled;  this  kind  of 
equation  is  called  the  SH  wave  equation  (Virieux  1984),  which  represents  a  veracity  in  the  plane  of 
symmetry  (mirror)  of  a  monoclinic  medium.  If  propagation  exists  in  this  plane  of  symmetry,  then  we  have 
anti-plane  strain  motion  which  represents  the  most  general  situation  under  which  pure  shear  waves  can 
exist  at  all  an^es  of  propagatioa  However,  hexagonal  media  possessing  pure  shear  wave  propagation  are 
deemed  to  be  part  of  a  degenerate  situation. 

If  we  consider  a  set  of  parallel  fractures  contained  integrally  within  a  transversely  isotropic  media,  then 
we  have  the  representation  of  a  monoclinic  medium;  here,  the  pure  anti-plane  strain  waves  are  SH  waves 
when  the  plane  of  symmetry  of  this  medium  is  vertical. 

In  addition,  it  can  be  said  that  other  cases  containing  symmetry  occur  in  monoclinic  media  such  as 
orthorhombic  media,  strong  trigonal  media,  and  weak  tetragonal  media 

If  we  consider  a  monoclinic  media,  Virieux  (1984)  suggests  that  the  elasticity  and  viscosity  tensors 
and  their  inverses  are: 
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ail  ^12  ^13  ®  hs  0 

^12  ^22  ^3  ®  ^5  ® 

^13  ^23  %3  ^  %5  ^ 

0  0  0  844  0  845 

^15  ®25  ^35  0  855  0 

0  0  0  845  0  Egg 


(25) 


Now,  according  to  Neumann’s  Principle  (Neumann  1885),  any  kind  of  symmetry  possessed  by  the 
attenuation  flows  in  conjunction  with  the  symmetry  of  the  crystallographic  network  of  the  substance  or 
material. 


Hence,  the  significant  elements  or  components  depicting  the  movement  of  the  SH  wave(s)  are: 


844  a4g 

^46  ^ 


(26) 


Thus,  it  can  be  readily  deduced  that  the  partial  differential  equations  can  be  extracted  from  the  second  row 
of  tensor  equation  (8)  and  the  fourth  and  sixth  rows  of  tensor  equation  (19)  as: 


3a^y  ^  80yg  =  o  -F 
d\  9z  dt  ^  ’ 


T44  Oy2  't4g  O^y  S44 


yz 


da. 


xy 


9t 


~^46- 


xy 


at 


(28) 


av  y  ^^XV 

=  '^46  ®yz  ®xy  +  S4g  ..  ...g.— +  Sgg  , 


(29) 
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where 


'^44  ~  ’  '^46  ®  ‘n46/^  » 


T1  =  1144^66  "  ^46  ’ 


(30) 


and 


S44  =  WC.  Sgg  =  C44/C, 


2 

c  -  C44  Cgg  -  , 


(31) 


such  that  the  stiffiiess  Cjj  and  the  viscosities  11^  (where  i,j  =  4,6)  are  the  (i*  j*)  compoiKnts  of  the  tensors 
Cm*  tIm’  respectively.  According  to  Carcione  and  Cavallini  (1993),  Equations  (22)-(24)  can  be  easily 
converted  into  Equations  (27)-<29)  reversibly  as  long  as  the  following  relations  are  taken  into  account 
(noting  that  s  and  t  can  be  introduced  as  simple  order  2  tensors): 


- 

- 

Hv 

y 

y 

^yz 

<=> 

-Ex 

1 _ 

(32) 


Fy  <=>  My  , 


S44  S46 

^11  -^13 

546  S66_ 

_-^13  ^3  _ 

T44 

'*^46 

Oil  -O12 

T  s 

<=> 

> 

'^66_ 

-<Ji3  O33  ^ 

P  ^  P> 


(33) 

e' ,  (34) 


o'  ,  (35) 


(36) 
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Similarly,  if  we  allow  the  stiffness  and  viscosity  tensors  to  be  of  order  2,  then  we  have, 


C46 

^66_ 


^44  ^46 

JI46  ^66 


(37) 


Hence,  we  can  arrive  at  the  2-D  identities 

S  =  C~^  and  T  =  Ti’^ , 


(38) 


typifying  the  3-D  Equations  (17)  and  (18),  respectively.  Therefore,  it  follows  logically  in  conclusion  that 
the  anisotropic  SH  wave  relation  premised  upon  a  Maxwell  theology  is  equivalent  to  the  anisotropic 
Maxwell  equations  in  a  mathematical  analogy  such  that  the  forcing  fimction  expression  would  be  the 
magnetic  current 
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Intentionally  left  blank. 
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